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ABSTRACT 



We report on the geometry of accretion disk and high energy coronae in 
the strong Comptonization state (the very high/steep power law/hard interme- 
diate state) based on a Suzaku observation of the famous Galactic black hole 
GX 339 — 4. These data were taken just before the peak of the 2006-2007 out- 
burst, and the average X-ray luminosity in the 0.7-200 keV band is estimated 
to be 2.9 x 10 38 erg s _1 for a distance of 8 kpc. We fit the spectrum with both 
simple (independent disk and corona) and sophisticated (energetically coupled 
disk and corona) models, but all fits imply that the underlying optically thick 
disk is truncated significantly before the innermost stable circular orbit around 
the black hole. We show this directly by a comparison with similarly broadband 
data from a disk dominated spectrum at almost the same luminosity observed 
by XMM-Newton and RXTE 3 days after the Suzaku observation. During the 
Suzaku observation, the QPO frequency changes from 4.3 Hz to 5.5 Hz, while 
the spectrum softens. The energetically coupled model gives a corresponding 
5 ± 8% decrease in derived inner radius of the disk. While this is not significant, 
it is consistent with the predicted change in QPO frequency from Lense-Thirring 
precession of the hot flow interior to the disk and/or a deformation mode of this 
flow, as a higher QPO frequency implies a smaller size scale for the corona. This 
is consistent with the truncated disk extending further inwards towards the black 
hole. 



Subject headings: accretion, accretion disks — black hole physics — stars:individual 
(GX 339 - 4)— X-rays:stars 
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Introduction 



The accretion flow in Black Hole B inaries (BHB) is generically unstable, so most 



sources are transient (e.g. 



Lasota 



200 ll ). showing dramatic outbursts with a large 



change in mass accretion rate and a correspondingly large change in X-ray properties 



(IMcClintock fc Remillard 



20061 ). Generally these start in the low/hard state (LHS) where 
the disk is very dim, and the spectrum is dominated by a hard Comptonized spectrum up to 
100-200 keV. This brightens rapidly while remaining hard, then there is an abrupt change, 
where the disk spectrum strengthens rapidly and the Compton tail softens to T > 2.0. The 
tail can then become quite weak, so that the spectrum is dominated by the disk emission 
with temperature < 1 keV of which the Wien tail extends into the RXTE bandpass above 
3 keV. This state is called the high/soft state (HSS). The source then remains in the 
HSS as its luminosity declines, then makes a transition back to the LHS. This typically 
happens at a lower luminosity than the LHS to HSS transition, a hysteresis effect which 
gives a charact eristic q shape on a hardness-intensity diagram of the outburst (see e.g., the 



compilation of 



Dunn et al. 



20101 ). 



In the HSS, t he disk emission can be fit by the diskbb model which approximates 



the standard disk (IShakura fc Sunvaey 



T(r) = T in • (r/r in )- 3 / 4 flMitsuda et al. 



19731) by describing the local disk temperature as 



19841 ) with the maximum observed disk temperature 



Ti n and an apparent disk inner radius r; n . In this modeling, the disk bolometric luminosity 
-^disk can be related to these two spectral parameters as L disk = 47irf n aT^. Compelling 
evidence for the standard disk formalism is given by the observation that the value of 
r; n remains remarkably constant in the HSS as Ldisk changes significantly. Thus, after 



several corrections inc 



1998 



Gierlihski et al. 



uding the stress-free inner boundary 



19991 ) , color temperature correction (IShimura &: Takahara 



and relativistic corrections (jCunningham 



1975 



condition (e.g., 



Kubota et al. 



Zhang. Cui fc Chen 



19951 ) 



19971 ) . r in is generally 



-4- 



believed to be consistent with the innermost stable circular orbit (ISCO) around central 
black holes. 

By contrast, in the LHS the disk emission is very weak and the luminosity is instead 
dominated by the hard (T < 2) power law tail. The radius of the disk is hard to measure 
directly in this state due t o model uncertainties as well as the weakness of the disk emission 



(see e.g 



Done et al 



20071 ). Models of this state generally assume that the inner disk 



progressively recedes at low luminosities, bein g replaced 



solutions of the accretion flow equations (e.g., 



EsinetaL 



jy th e alternative hot inner flow 



19971 ) . These models then predict 



that as the source makes a transition from the HSS towards the LHS, the disk should start 
to pull back from the ISCO. This cannot be observed directly in RXTE data due to the 
3 keV lower limit to the bandpass. However, there is one source monitored by SWIFT 
during a transition, and here the radius of the thermal disk component cle arly increases as 



predi cted as the the source declines from a HSS into the intermedia te state ( 



200 81), although the LHS radius is much more model dependent (IRykoff et al 



Gierlihski et al. 



2008|). 



Gier 



ihski et al. 



2007 



The intermediate states from LHS to HSS on the rise often appear quite similar to the 



spectra seen as the source makes a transition f rom HSS to L 



lumin osity due to the hysteresis noted above (iBelloni et al 



IS on the decline, but at highe r 



1996 



Mendez van der Klis 



19971 ). Some of the intermediate spectra have st eep power law tails with r > 2 .4 leading 



to them being called the steep power-law state ( IMcClintock fc Remillard 



20061 ). although 



a power law is not a good approximation for the complex curva ture of Comptonization 



seen in this state ( Zdziarski et al. 



2002 



Gierlihski fc Done 



20031 ). This state is much more 



common on the rise than the decline, so is generall y a high luminosity stat e, leading to its 



original name of a very high state (hereafter VHS, 



Miyamoto et al. 



199 lh . Where these 



spectra are seen on the decline, the disk is always clearly seen as a separate component from 
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the Compton tail. However, some of the most extreme examples of this state seen on the 
rapid rise are simply dominated by the Comptonized component, sometimes without even 
a point of inflection in the spectrum to mark the presence of the disk (e.g. RXTE data of 
ObsID 30191-01-02-00 of XTE J1550 - 564 and 91702-01-58-00 for GRO J1655 - 40). 

The different states can also be identified by their fast variability properties. In the 
HSS the disk spectrum is remarkably constant on timescales less than a few hundred 
seconds. This is expected as the viscous timescale is long for a geometrically thin, cool disk. 
Thus there is very little variability power at frequencies / < 0.01 Hz at energies where the 
disk dominates the spectrum. Conversely, in the LHS there is strong band limited noise, 
with a quasi-periodic oscillation (QPO) superimposed. The frequency of this QPO termed a 
type C increases from 0.1-6 Hz as the spectrum softens from the LHS into the intermediate 
states. These intermediate states can themselves be split into a hard intermediate state 
(HIMS) and soft intermediate state (SIMS). As suggested by the nomenclature, the SIMS 
has stronger disk and weaker tail than the HIMS, but this is a rather subtle distinction. 
The main difference between these two states is that the HIMS has variability properties 
which are clearly similar to those in the LHS, with broadband noise and type C QPO, 
whereas the broadband noise collapses in the SIMS, leaving the power spectrum dominated 
by the QPO alone (termed a type B QPO where this still retains the harmonic structure 



seen in the HIA. 


S, or 


Wiinands et al. 


1999: 



Belloni et al. 



2002 



Belloni 



20101 ). The most extreme VHS spectra 



noted above have power spectra like the HIMS, while those where the disk is clearly visible 
as a separate component most often have power spectra which are similar to the SIMS. 

The high energy emission in the VHS, HSS, and bright LHS, is caused by energetic 
electrons via inverse Compton scattering of low energy seed photons from the disk. In the 
HSS, this corona is either optically thin or patchy, so only intercepts a small fraction of the 
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disk flux, making very little difference to the observed r; n . However, in the VHS/HIMS 



and VHS/SIMS, t 



re corona is clearly 



of the inner disk (IZdziarski et al. 



3 0th moderately optically thick and covers mos t 



2002 



Gierlihski fc Done 



2003 



Kubota fc Done 



20041 ). 



Compton scattering removes photons from the disk, boosting them in energy to form the 
tail. Hence the disk appears less luminous, it has an apparently smaller radius (Kubota 
et al 2001; Steiner et al 2010). The constant HSS radii can be recovered by correcting 



for this effect in spectra where the corona carries less than 



( IKubota. Makishima fc Ebisawa 



2001 



Kubota fc Makishima 



lalf of the accretion pow er 



2004 



Steiner et al. 



20101 ). but 



beyond this point the inferred radius of the disk is larger than in the HSS, indicating that 



the disk may be trun cated before the ISCO (IKubota fc Done 



2004 



Done fc Kubota 



2006 



Gierlihski et al. 



20081 ). It is essential to understand the precise geometry of the optically 
thick disk and X-ray coronae in these states in order to derive constraints on the inner disk 
radius and hence test whether the disk begins to recede as predicted by the alternative hot 
accretion flow models. 

In this paper, we utilize Suzaku data of a well known stellar black hole GX 339 — 4 in 
the VHS/HIMS seen during the 2006/2007 outburst which is the brightest outburst ever 
observed by RXTE. In the energy range of 1.5-12 keV, the RXTE/ASM count rate of 
peak of the 2006/2007 outburst was ~ 14% higher than that of the 2nd brightest outburst 
in 2002-2003. This is one of the best available VHS/HIMS spectra, and hence the best 
to constrain the disk geometry during the VHS/HIMS in a LHS-HSS transition. Several 
parameters of GX 339 — 4 and its outburst during 2006-2007 are summarized in §2. The 
observation and data reduction are briefly described in §3. The Suzaku spectra are analyzed 
using the commonly used diskbb plus power-law model (§4.1), two independent corona 
models (§4.2), and an inner disk-corona coupled model (§4.3), and then, the truncated disk 
is discussed in §5. In §6, QPO frequency and change of the disk parameters are briefly 
discussed. Errors quoted in this paper represent 90% confidence limit for a single parameter 
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unless otherwise specified. 



2. GX 339 — 4 and its X-ray properties 



GX 339—4 is a transient Galactic black hole, discovered in the early 1970s flMarkert et a 



973). The system parameters are rather poorly known from optical studies. 



Hynes et al. 



(120041 ) give a lower limit of 6 kpc, but their favored distance of 16 kpc is not consistent 



with the measured ab sorption systems in GX 339 — 4, so a more likely distance is 7-9 kpc 



(IZdziarski et al. 



20041 ). This matches well with models of the binary system, where the 
low mass companion star must be big enough to have a high mass transfer rate to power 
the observed lumin osity so needs to be substa ntially more distant than 6 kpc so as not 



to be seen directly (jMunoz-Darias et al. 



20081 ). Hereafter we assume a distance of 8 kpc. 



The orbital inclination is estimate d to be less than 75° to the line of sight, as inf erred 



from the lack of eclipses and dips (jCowley et al. 



2002 



lack of strong iron absorption lines in the HSS also argues for % < 60° ( 



Kolehmainen & '. 


-I 

Done 


2010). 


The 


ues for i < 60° f 


Ponti et al. 


2012) 



Conversely there are lower limits on the inclination from studies of the binary (i > 45°: 



Hynes et al 



so argues for ■ i > 50° if this 



20031 ) while the strong QPO seen in this source a 
is made from Lense-Thirring (vertical relativistic) precession (jlngram et al.ll2009l ). We use 
i = 50° as it fills all these constra ints and is the result obtained from the iron line profile 



analysis of IShidatsu et al. 



(2011b) 



The HSS spectrum of this source was for the first time obtained with Tenma, and 



Makishima et al. 



(119861 ) estimated the apparent inner radius. With our system parameters 
as used above this corresponds to r in ~ 52 ■ Cso^s km, where C50 is (cosi/ cos 50 )- 1 / 2 with 
the inclination angle i, and dn is th e distan ce to the source in the unit of 8 kpc. Similar 



results are found by 



Shidatsu et al. 



(2011a 



based on the 2-9 keV MAXI and the 0.6-7 



keV Swift data during the 2010 outburst. The RXTE study of 



Kolehmainen fc Done 
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( 120 lOl ) with multiple datasets gives a larger estimate, r in as 67.5 ± 0.9 • Cso^s km. 
Changing the bandpass to the lower energy range of CCD data gives a slightly different 
normaliz ation as the diskbb mod el does not completely describe the spectral shape of 



the disk ( IKolehmainen et al. 
give r in = 57 • Cso^s km (iKolehmainen et al. 



2011). The three disk dominated spectra in XMM-Newton 



201ll ). Thus the value of apparent inner 



radius in the HSS is estimated as r in ~ (52 ~ 68) • Cso^s km. Folding in correction factors 
of color temperature, k, and inner boundary condition, £, gives a radius for the ISCQ 
of -Rt sco = rin • ^ 2 • £ — (56 ~ 73) • Cnncfa km, taking k ~ 1.7 (jShimura &: Takahara 



19951 ) and £ = 0.37 (IGierlihski et al. 



19991 ). The alternative approach of fitting proper 



radiative transfer disk mode ls with full relativistic corrections also gives similar values 



(IKolehmainen fc Done 



20101 ) . 



Shidatsu et al. 



(I2011bh . 



The disk structure in the LHS has been analyzed in detail by 
They used Suzaku data of this source obtained in 2009 March, and showed that the iron line 
features indicated the inner radius of 13.3^6 Rgi where R g = GM/c 2 is the gravitational 
radius. It corresponds to 140^ (M/7M ) km. They also showed that weak continuum 
emission from the optically thick disk and surrounding coronae was consistent with the 
truncated disk picture with Ri n ~ 120 • Cso^s km, though again we caution that this is much 
more model dependent than the HSS. 



3. Observation and Data reduction 



As described in §1, the 2006/2 007 outburst of GX 339 — 4 was its brightest outburst 



since the launch of RXTE in 1995 ( jSwank et al. 



20061 ). Figure ffl shows the 1.5-12 keV and 



15-50 keV light curve and hardness ratio of this outburst, obtained with the RXTE/ASM 
and the Swift /BAT Transient Monitor. The time histories of the hardness ratio and the 
Swift/BAT count rate suggests that the state transition occurred around MJD~54140- 
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54150. The Suzaku observation was performed on 2007 February 12th 05:33:31 to February 
15th 04:48:26 corresponding to MJD=54143. 2-54146. 2 as indicated with dashed lines in 
figure [TJ Clearly the source hardness and the BAT count rate changed dramatically around 
the time of the Suzaku observation. An analysis of the QP O seen in RXTE showed a change 



from type-C to type-B on February 16 ( IMotta et al 



20091 ). so the source made a transition 



from VHS/HIMS to VHS/SIMS between February 15 and 16. The Suzaku spectrum and 
variability shows this to be a VHS/HIMS spectr um, though ; 



marking the separation of disk and Compton tail (IMiller et al 



2008 



3 still an inflec 


tion 


Yamada et al. 


2009) 



The XIS employed the 1/4 window option and a 
and 0.5 s for XIS3). The Hard X-ray Detector (HXD; 



burst option (0 3 s fo r XIS0/XIS1, 



Kokubun et al. 



20071 ) was operated 



in the standard mode . The data processing and reduction were performed in the same 



way as 



Yamada et al. 



(120091 ) . using the Suzaku pipeline processing version 2.0.6.13. To 
extract light curves and spectra, the data were analyzed with HEAsoft version 6.9 and the 
calibration data files (CALDB) released on 2007 July 10th. The XIS events were extracted 
from a circular region with a radius of 7' centered on image peak. Since this extraction 



circle is larger than the window size of 17'. 8 x 4'. 5, the effective extraction 
the in tersection of the window and this circle. As reported in detail by 



region is there 



Yamada et al. 



ore 



( 120091 ). the XIS events piled up significantly, and thus a central region of 3' at the image 



center was exc l uded from the event extraction region to minimize these effects. Following 



Yamada et al. 



(120091 ). we use data from XIS0, since data from XIS1 and XIS3 are affected 
more by pile up than XIS0. In addition, a variable fraction of the CCD frame was often 
lost due to tele metry saturation , resu lting in only 2.83 ks of data out of the 12.4 ks of XIS0 



exposure. While 



Yamada et al. 



( 120091 ) included data affected by the telemetry saturation to 
describe the iron line in detail, we excluded the telemetry saturation to avoid a small effect 
on the continuum shape due to the asymmetric response. 
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The PIN and GSO spectra, acquired for a net exposure of 87.9 ks, were corrected 
for small dead time, but no other correction due to the source br ightness was 



neces sary. We subtracted modeled non-X-ray backgrounds (NXBs; 



Fukazawa et al. 



20091 ). The PIN and GSO background spectra were constructed based on lcfitdt 
and lcfit methods (2.0ver0804), respectively, filtered by the same good time intervals. 
As the response files, ae_hxd_pinxinome3_20080129.RSP was used for PIN, and 
ae_hxd_gsoxinom_20080129.rsp and ae_hxd_gsoxinom_crab_20070502.arf were 
used for GSO. The cosmic X-ray background was ignored, since it is less than 1% of the 
total counts. We use GSO data only up to 200 keV, since the signal becomes smaller than 
3% of the background at 200 keV. Figure [2] shows the Suzaku light curves of GX 339 — 4 in 
which background levels were subtracted and dead times were corrected for PIN and GSO 
data. 



4. Analyses of averaged spectrum 

In this section, we present the analysis of the Suzaku spectra extracted above. To fit 
the observed spectra, we use the xspec spectral fitting package (version 12.6.0). We use 
energy range 0.7-9.0 keV for XIS0, 13-60 keV for the PIN spectrum, and 70-200 keV for the 
GSO spectrum. Large fit residuals due to calibration uncertainties are often observed near 
the edge structures of the XIS/XRT instrumental responses so we exclude XIS data from 
1.4-2.3 keV. We fit the three spectra simultaneously, with constant factors scaling between 
the three fixed at 1, 1.07, and 1.07 for XIS, PIN and GSO, respectively^. We extend the 
energy range of the spectral fitting to 0.1-1000 keV as we use some convolution models. 
Such models have edge effects at the end of the energy range used for the calculation, so it 



^ttp:/ /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/watchout. html 
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is important that this beyond the energy range used for the data. 



4.1. Empirical modeling by DISKBB and POWER-LAW with reflection (model 1) 



In order to characterize the spectral shape , we fi rst model the 0.7-200 keV 



Suzaku data in the same way as lYamada et al 



(120091 ). i.e., the diskbb model plus 



power-law (hereafter PL) model. We hereafter call this 'model 1'. The two continuum 
components are both absorbed by a c ommon hydrogen column mod eled with xspec wabs 



model (IMorrison fc McCammon 



Yamada et al. 



(120091 ) used LAOR model to 



1983|). While 

describe the iron line in detail, we use a single gaussian with a at 0.2 keV since we focus on 
the continuum shape and our much shorter exposure (due to more stringent constraints on 
telemetry issues) means that the spectra have much lower statistics. 

To account for ionized reflection of the PL compo nent, we utilized the model 



ireflect (a convolution version of the pexriv model 



Magdziarz fc Zdziarski 



19951 ) 



This model balances photoi onisation with ra diative recombination from the very simplistic 



ionization balance code of (IDone et al. 



19921 ). This requires a user defined (rather than 
self-consistently computed) temperature, which we set to 10 6 K as the disk is clearly hot. 
However, the code does not calculate collisional processes, so the ion population are set 



only by photoionisation. This is clearly a p oor assumption 



fairly high ionisation parameter required by 



Yamada et al 



or such a hot disk. Thus the 



(120091 ) may be an overestimate 



of the illuminating flux required. Yet another limitation is that Compton scattering within 
the disk is not included, which has quite a large effect on the characteristic iron line and 
edge profile in the reflected emission. Photo-ionization heats the disk photosphere to some 
fraction of the Compton temperature, so this can be of order ~ 1 keV for high ionization 
parameters. All the line/edge features are produced at an optical depth of around unity, 
so this means that a third of the photons are Compton scattered by these hot electrons, 
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smearing the sharp line/edge seen in ireflect/pexriv (iRoss et al.lll999l ; iRoss fc Fabian 



20051 ). This may mean that the relativistic smearing is overestimated, but we caution that 
generally spectral parameters are all interrelated, so it may instead distort the observed 
solid angle, ionization state or iron abundance. Despite these drawbacks, we use this as 



there are few better alternatives at present. The models at high densities of 
(120071 ) are not yet public. 



Ross fc Fabian 



We fix the ionization parameter 300 erg cm s _1 , and fix all abundances other than 
iron at solar. Hence the only free parameters are the amount of solid angle and the 
iron abun dance. We smear this reflected spectrum by general relativistic effects using 

. We fix all the parameters of this at characteristic values 



rdblu r (|Fabian et al. 
seen by lYamada et al 



1989 



f]2009h . i.e. inner and outer radius R^ bluv = !0R g , R^ lur = 10 5 R 9 , 
inclination angle % = 50°, and power law index of emissivity = —3 as our limited data 
quality do not allow us to constrain them independently. 

Table [1] gives details how the model was described in xspec, while Table [2] gives the 
best fit parameters. This shows that model 1 fits the data well with x 2 /dof = 190.2/186. 
Figure [3] shows the data and the best fit model and residuals, while figure 0] and figure [5^ a) 
show the vF v spectrum and unabsorbed model spectrum, respectively. The absorbed 
0.7-200 keV flux is estimated as 2.8 x 10 -8 erg s _1 cm -2 which gives an absorbed luminosity 
of 2.1 x 10 38 d% erg s _1 assuming isotropic emission. The PL photon index Tpl = 2.68lnm , so 



this p lus the high luminosity means that the source is in the VHS, as n oted by 



Yamada et al. 



Yamada et al 



(1200 



( 120091 ). The best fit parameters are consistent with those obtained by 
within 90% confidence, and the best fit diskbb model shows disk inner temperature of 
kT in = 0.65 ± 0.02 keV and an apparent inner radius of r in ~ (56^g) • Cso^s km, consistent 
with that observed in the usual HSS, (52 ~ 68) ■ Cso^s km (see §2). 



However, an inspection of the figure [5|a) shows that this is not a good physical 
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description of a model where the disk provides the seed photons for Compton up-scattering 
into the PL tail as the PL extends below the disk at low energies. This motivates further 
studies of the disk structure with more physical models. 



4.2. Disk geometry based on independent corona (model 2 and model 3) 



Compton scattering conserves photon number, so it removes seed photons from the 
observed spectrum by boosting them in energy to form the Comp ton tail. We replaced 



20091 ) which self consistently 



the PL model with a convolution model simpl (jSteiner et al. 
removes as many seed photons as it puts into a Compton tail. This is identified as 'model 2' 
in table [TJ As described in table [TJ only the up-scattered photons are reflected by the 
optically thick disk, and thus modified by the models of rdblur and ireflect. The fit 
results are presented in table [2] and figure [3](c). As shown in unabsorbed model spectrum 
(figure Eb), this gives a more physical description of the distribution of the low energy 
photons. The lack of emission at the lowest energies from Comptonization compared 
to a PL means that the disk temperature becomes lower to compensate for this, with 
kT in = 0.551q;o^ keV, and the radius becomes higher r in = (122+g) ■ Cso^s km. This implies 
that the disk inner radius is twice as large as the ISCO, by comparison to that found in the 
usual HSS. 



However, the Comptonization se en in the VHS is 



mal and non-thermal electrons (e.g. 



comp 



Gierlihski et al. 



1999 



ex, containing a mix of both ther- 



2001 



Gierlihski fc Done 



Kubota. Makishima fc Ebisawa 



20031 ) whereas the simpl model used above only makes a 



non-thermal power law tail. Hence we a dded a thermal comptonization model, nthcomp 



( iZdziarski et al. 



1996 



Zycki et al. 



19991 ) to the model (model 3 in table [Q. This has 
four parameters, seed disk photon temperature, which we tie to that of the seed photon 
temperature of diskbb, the electron temperature kT e , photon index r t h (equivalent 
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to optical depth), and normalization. The relation betwee n the optical depth, e lectron 



Zdziarski ct al. 



(Il996h assuming 



temperature and spectral index is given by equation (Al) of ! 
a spherical source with a uniform distribution of seed photons throughout the source. By 
contrast, we assume slab geometry with seed photons at the botto m of the slab. Comparin g 



Poutanen Svensson 



these two geometries with compps (geom=— 4 and 1 respectively, 
19961 ) shows that the spectral indices are equivalent for an optical depth in the slab 
geometry which is approximately half that of the sphere. Hence we use 



1 

T= 2 



- 



4 e c • ((r th + \y - 



3 



(1) 



where G e = kT e /m e c 2 . 



We cannot uniquely constrain both T t h from the thermal Comptonization, and T p i 
from simpl so we fix the la tter at 2.1, the value of the non-thermal tail seen in the typical 

19991 ). T he non-thermal tail may be somewhat steeper in the VHS, 



HSS (jGierririski et al 



and again has a complex shape (IGierlihski fc Done 



20031 ) . However, our data do not extend 



above 200 keV so we are not very sensitive to this. Hence model 3 has only two additional 
free parameters. Similarly to models 1 and model 2, only the Comptonized emission is 
reflected by the underlying optically thick disk, and thus the model description is somewhat 
complicated as described in tabled) This gives y^jdof = 192.0/184, still not as good as 
model 1, but better than model 2. The fits are shown in tabled figure [3]^d) and figure [5](c). 
These show that the thermal Comptonization dominates the spectrum below ~ 50 keV. 

Unlike simpl, the thermal Comptonization does not decrease the seed photon 
spectrum. Hence we have to manually adjust the disk normalization for the photons 
which are scattered out into the thermal Compton tail. Hence we also tabulate the 
unabsorbed photon flux jp^ oton contained in the thermal Compton tail integrated from 
0.01-100 keV. We assume that these photons were removed equally from all energies of the 
disk spectrum (equivalent to a uniform corona over the entire disk). The observed photon 
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flux from the disk, i^£ ton is a l so tabulated, and we simply increase the measured disk 
normalization by this amount. Since the normalization is proportional to the square root 
of the radius this gives an inferred radius r* n = r in ^J (^F^ oton _|_ F^f on )/Fj*£ ton , where r in 
is calculated from the diskb b normalization alone. This is equivalent to equation (A.l) in 
Kubota fc Makishimal (120041 ) . and this gives r* n = (117 ± 5) ■ Cso^s km, which is as large as 



that estimated with model 2. 



4.3. Inner disk-corona: coupled energetics (model 4) 



The previous two models assumed that the Comptonizing electrons fully covered the 
disk, and that the underlying disk structure is unaffected by the presence of the high energy 
electrons. However, both these assumptions are probably not appropriate. The corona 
probably only covers the inner disk, at r in < r < r tran as schematically shown in figure [6j 
The power in the corona must also derive ultimately from the accretion flow, so there 
should be less power dissipated in the disk. 

In the view of these concepts, we fit the data by replacing the 'diskb b+nthcomp' 



Done &: Kubota 



compo nents in model 3 to the coupled disk-corona model, dkbbfth given by 
(120061 ) (model 4 in table [1]). This model assumes that the energy released by gravity is 
dissipated locally, either thermalising in the disk for r > r tmn , or split between the disk and 
corona for r m < r < r tran , as shown schematically in figure EJ Thus the disk underlying the 
corona is cooler and less luminous than it would have been if the corona did not exist, and 
it is this weak a nd cool disk emission which provides the seed pho tons for the Compton 



upscattering (see 



Svensson fc Zdziarski 



1994 



Done Kubota 



20061 ). 



The model parameters are similar to that of using a diskbb+nthcomp continuum 
(model 3), except that rather than having two separate normalizations, the parameters 
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which control the relative normalization of the disk and corona are a combination of 
r tram which determines the outer, unComptonized disk luminosity, and the shape of 
the thermal Comptonization i.e. r t h (equivalent to r) and kT t ^. The luminosity of the 
Compton component, L t h, is approximately determined by the Compton y parameter 
where y « 46 e (r + r 2 ) as L th = yL disk - m and 

-^diskjin is the inner disk luminosity for 
r in < r < r tran . Thus the ratio of power dissipate d in the corona to tha t of the disk in the 



region r in < r < r tran is / th = L t h/L diskiin « y (see 



Done fc Kubota 



2006 



for examples of how 



the model spectra cha nge as a function of th e se pa rameters) . The inner radius is calculated 



Kubota k Makishimal (120041 ) by replacing 'F d p isk + F t p hc • 2 cos V to 



via equation (A.l) in 

photon flux of the dkbbfth component, i^dkbbfth; assuming the slab geometry. The radius is 
also determined from the intrinsic (rather than observed) innermost disk temperature, T; 1 "*. 
This is how the entire disk would be seen if there was no energy dissipated in the corona. 

This model only takes into account the thermal Comptonization, so we again include 
the simpl with fixed Tpl = 2.1 to describe the non-thermal tail. In this analysis, we set the 
un-scattered component of the dkbbfth to be seed photons for simpl. Again, both the 
thermal (dkbbfth) and non-thermal (simpl) Compton components are both modified by 
the ionized reflection and blurred by the relativistic effects. 

The model has the same number of free parameters as model 3, but gives an improved 
fit with x 2 /dof = 178.7/184. This is the best chi-squared value among the four models 
despite the fact that the model is more constrained. The fit results are summarized 
in table [21 figure [3](e) and figure El^d). With the best fit parameters by assuming an 
isotropic emission, the absorbed and unabsorbed luminosity was estimated in the range 
of 0.7-200 keV as (2.1 x 10 38 )d 8 erg s" 1 and (2.9 x 10 38 )d 8 erg s" 1 , respectively. We 
also calculated the bolometric luminosity (3.8 x lO 38 )^ erg s _1 and the optical depth of 
r = 0. 58^g'o 8 via equation (1). The corona was found to be localized within r tran of 2.7 r* n , 
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within which, 27% of the accretion power is dissipated in the thermal corona. 

Considering the power dissipated in the corona, the intrinsic disk temperature should 
be kTg* = 0.67±o;o8 keV - With this intrinsic disk temperature and the 0.1-100 keV 
DKBBFTH photon flux of ^dkbbfth = ^3.5 photons s _1 cm -2 , the underlying disk inner radius 
was estimated as r* n = (93tg 9 ) • Cso^s km. In figure 0(d), the intrinsic disk emission is 
plotted together with a dash-dot line. Though the estimated apparent inner radius is 
slightly smaller than that obtained under the independent corona modeling (model 2 and 
3), it is ~ 1.3-2.2 times as large as that observed in the HSS (see §2). 

We checked that the results were not affected by the assumption that relativistic 
blurring has a radial dependence parameterized by /3 = —3. We replaced this by (3 = 10, 
which is more appropriate for a stress-free inner boundary condition by assuming that the 
line emissivity is proportinal to the local energy release rate in the disk. The difference 
from the (3 = —3 case was negligible. Even removing relativistic effects entirely does not 
change the continuum parameters significantly due to the limited statistics around the iron 
line/edge in our data. In all cases the best fit disk parameters were still consistent with 
those shown in table [2j 



Is the disk truncated in these data? 



All the models above which self consistently correct the disk normalization for Compton 
scattering (models 2, 3 and 4) find a disk inner radius which is larger than that seen in 
the HSS in this source. This can be seen directly by comparison with a HSS spectrum 
at a very similar luminosity. Figure [7(a) shows the Suzaku data used here together with 
the H SS spectrum of almost th e same luminosity taken by XMM and RXTE on February 



19th (jKolehmainen et al. 



201 ll ) (see figured]). Figure UJb) shows the comparison of the 
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unabsorbed model spectra. 

The HSS spectrum is roughly characterized by a dominant diskbb of kT in ~ 0.82 keV 
with a weak PL tail, and an absorbed 0.7-200 keV flux is estimated as 2.4 x 10~ 8 erg s -1 cm™ 2 , 
which is only 15% lower than the Suzaku flux of 2.8 x 10~ 8 erg s _1 cm -2 in the same energy 
band. Thus if the disk has constant radius between the HSS and VHS/HIMS then the disk 
temperature should be slightly higher in the higher luminosity VHS/HIMS. Yet is it clear 
from this figure that the VHS disk has a lower temperature, as the data peak at lower 
energy in figure [7(a). Compton scattering retains the imprint of the seed photon energy, 
so this lower temperature c annot be a consequenc e of simply Comptonizing the higher 



temperature inner disk (see iDone fc Kubotal 120061 ). The blue dash-dot line in figure UJb) 



shows the estimated intrinsic disk emission from the VHS/HIMS data as reconstructed from 
the dkbbfth model, i.e. how the disk emission would have looked in the Suzaku data, if 
the thermal corona was not present and the matters of thermal corona were all accreted in 
the disk. If the inner radius is kept constant between the HSS and the VHS/HIMS, the 
intrinsic disk emission should peak at higher energy than in the HSS. However, the peak 
energy of the intrinsic disk emission is still much lower than that in the HSS. Estimated 
intrinsic VHS disk temperature kT^ = 0.671qo8 keV * s mucn lower than the observed 
HSS disk temperature of kT[ n ~ 0.82 keV. Therefore, the apparent inner radius of the 
VHS/HIMS is clearly larger than that in the HSS. 

To convert from apparent to true inner disk radius requires corrections for the color 
temperature, k, and inner boundary condition, £, with R in = r in ■ k 2 ■ £ as described in §2. 
If we use the same n and £, we obtained R m = (99lg°) ■ Cso^s km. Under the same value of 
£, the only way to get the same inner radius of the disk in the VHS/HIMS as in the HSS is 
if k is as low as ~ 1.3-1.4. This seems very unlikely given the strong illumination. On the 
contrary, it is certainly plausible that the disk color temperature correction has increased 
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substantially, as the inner disk should be strongly illuminated by the Comptonized emission 
in the VHS/HIMS whereas it is not in the HSS. This only reinforces the change in true 
inner radius as a larger color temperature correction means that -R; n is even larger than 
derived from assuming k = 1.7. Similarly, if the disk in truncated, then the use of the same 
boundary condition is not appropriate as it now does not extend down to the ISCO, where 
there is the stress-free inner boundary, but is truncated where there can still be stress on 
the inner bou ndary, making £ large r. This again leads to an increase in i?i n (see also the 



discussion in (jGierliriski et al. 



20081 ). 



Thus both of the expected changes in color temperature and inner boundary condition 
reinforce the conclusi on that the disk inner radi us is larger in the VHS/HIMS than in the 



HSS. Using k = 1.9 (IShimura fc Takahara 



19951 ) and £ = 1 as a reference, R in can be as 



large as ~ 340 • Cso^s km. 



6. Variability of the disk and corona 

6.1. QPO frequency 

Figure [8]^a) shows the power spectral density from the entire PIN lightcurve calculated 
using xronos (version5.22). In this figure, a double QPO feature seen at frequencies of 
~ 4 Hz and ~ 5-6 Hz. During the observation, PIN and GSO count rates slightly decreased 
by ~ 25% and ~ 18%, respectively, while that of XIS-0 slightly increased ~ 9% (figure [2]). 
Thus the double QPO feature is probably due to a change of a single QPO frequency. 

We reanalyzed the same data by dividing them into first half and latter half as indicated 
in figure |2j The calculated power spectral densities were also shown in figure E^b)(c) showing 
a clear change in QPO frequency from 4.26 ± 0.06 Hz in the first half to 5.46 ± 0.06 Hz, in 
the second half. Since these QPO frequencies are probably set by the size of the emission 



-20 - 



region, we look for changes in the size of the disk and corona in their corresponding 
spectra. Specifically, a higher frequency QPO should indicate a smaller size region, and our 
modeling with dkbbfth means that we can track both the size of the disk and the size of 
the corona via the parameters r* n and r tra n- We then compare the observed changes with 
the predictions of two specific QPO models. 

6.2. Time evolution of spectral parameters 

To investigate evolution of the disk and corona in detail, we accumulated spectra from 
the two halves of the observation. For PIN and GSO data, net exposures of 43.9 ks and 
44.0 ks were acquired for the first half and the latter half, respectively. The NXB spectra 
were calculated for each periods and subtracted from the data. For XISO data, 1.6 ks and 
1.3 ks exposures were acquired for the first half and the latter half, respectively. The two 
data sets have very similar observed luminosities. The absorbed flux was almost same as 
2.7 x 10 -8 erg s _1 cm -2 and 2.8 x 10~ 8 erg s _1 cm -2 , for the first half and the latter half, 
respectively, while their spectral shape was clearly different. The top panel of figure [9] shows 
vF v spectra of the first half (green) and the latter half (orange), and the bottom panel 
shows ratios of each spectrum to the best fit model 4 for the summed spectrum. There is a 
clear anti-correlation between the soft band data and the hard band data. The XISO data 
in 1-4 keV increased, while the PIN data decreased, so the spectrum pivots around 10 keV, 
becoming significantly softer during the Suzaku observation. 

A closer inspection of figure reveals that the data below 1 keV do not take part in 
this spectral pivoting, but remain relatively stable, despite the 10-20 % flux change in the 
range of 2-4 keV. These photons below ~ 1 keV are from the outer part of the disk, where 
the local temperature is lower than ~ 0.3 keV. Since this emission does not change it seems 
most likely that the mass accretion rate through the disk is not changing (as is also implied 



-21 - 



by the constant luminosity), and it is only the geometry of the inner disk/corona which 
gives the spectral change above 1 keV. 

In order to discuss change of the spectral parameters, we fit the two data sets with the 
independent corona (model 3) and the coupled corona (model 4). We fixed the values of 
iVn, iron abundance, and the gaussian central energy at those of the best fit values obtained 
by the summed spectral fit (table [2]). Furthermore, in the case of model 4, to avoid a strong 
coupling between r and kT e , we fix kT e to the value seen in the average spectrum. 

The results are summarized in table [3J and again the energetically coupled corona 
model dkbbfth gives the best fit. Based on the best fit model 4, the photon index of 
thermal corona r t h increases from 2.59 ± 0.03 to 2.74 ± 0.04. Since we keep kT e constant 
between the two datasets, this means that the optical depth of the thermal corona r slightly 
decreased from 0.61 ± 0.02 to 0.55^q q^. Hence the second half of the data have a lower 
y parameter and softer Compton tail, as seen in Figure |9^a). This means that / t h also 
decreases, from 0.30 to 0.25. The unComptonized disk emission is slightly more prominent 
in the second spectrum, so the transition radius decreases slightly from r tran = (2.6^2)?"^ 
to (2.3l^)r*, 

The intrinsic disk temperature kT^ marginally increases from 0.66 ± 0.03 keV to 
0.69+ 

ao3 keV, giving an apparent inner radius of the underlying disk, r * n which decreases 
by ~ 5% from (95+g) • Cso^s km to (90^) ■ Cso^s km. Thus the outer radius of the corona, 
^tran decreases by ~ 16% from 250 • Cso^s km to 210 • Cso^s km. Model 3 gave a similar result 
that the disk inner radius decreased ~ 5% from the first half to the second half data. 
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6.3. Discussion on the relation of QPO and the thermal corona 



Lense-Thirring preces sion has been proposed as the origin of the low-frequency QPOs 



(e.g. 



Stella fe Vietri 



1998). This can explain both the frequency and the spectrum of 



the QPO if the prec ession is of a hot inner Comptonizing flow rather than the thin disk 



( llngram et al. 



20091 ) . In this model the QPO frequency /qpo is /qpo r 21 . Thus, a 



28 ± 3% increase in QPO frequency means that the size of the precessing region should 
decrease by ~ 11 ± 1%. Since the precession is a vertical mode, the only part of the corona 
which can precess in this way is the corona interior to r* n (see the schematic diagram in 
figure [6]). The fits imply that r* n decreases by ~ 5 ± 8%, consistent with the Lense-Thirring 
QPO model. 

Another model for the low frequency QPO is that it is a mode of the hot inner flow. 
Global 3D magnetohydrodynamic (MHD) simulations suggest that the magnetorotational 
instability (MRI) m ay lead to the quasi-periodic deformation of a hot inner torus from a 



circle to a crescent (IMachida fc Matsumoto 



20081 ). The frequency in their simulation can 



be roughly related to the size of the inner torus as /qpo oc r~ 3 / 2 . So the observed change 
of QPO frequency requires that the radius change by ~ 15%. This is slightly larger than 
the observed decrease in r* n . However, their simulation was for an unconstrained hot flow, 
i.e. there was no cool disk present. It is not clear how this will change the formation and 
deformation of such a hot inner torus. So the model is probably consistent within the 
systematic uncertainty. 



7. Summary and conclusions 



We have analyzed the Suzaku spectrum of GX 339 — 4 in the strongly Comptonized 
VHS/HIMS. We used a series of models to fit the data, starting with the commonly used 
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'diskbb plus power- law' model (model 1). We then use progressively better models of 
the Comptonization, first assuming that it is non-thermal (model 2) and then allowing 



it to be both t 


lermal and non-thermal 


mode 


Zdziarski et al. 


2001 




Gierlihski & Done 


2003) 



20031 ) . The strong Comptonization means that 



many seed photons from the disk are removed by Compton scattering so the disk luminosity 
is larger than the directly observed luminosity. This effect increases the apparent radius 
of the disk, though using proper Comptonization models rather than a power law already 
increases the disk luminosity and gives an apparent radius which is larger than that seen in 
the HSS. Our final model couples the disk and corona together (model 4). 

Both these models except for model 4 assume that the disk and corona are independent 
components, yet both ultimately derive from the gravitational energy. Hence our final 
model is a more physical one where an inner corona reduces the power available for the 
inner disk. This more physical model gives the best fit to the data, and again requires 
that the apparent radius of the disk is larger than that seen in the HSS. This increase in 
apparent size of the disk can be seen directly by comparison with an XMM-Newton/RXTE 
HSS spectrum at almost the same luminosity taken 3 days after the Suzaku data. This 
clearly shows that the disk temperature is lower in the VHS/HIMS than in the HSS. The 
only way that this can be consistent with the same true inner disk radius is if the color 
temperature correction and/or stress at the inner boundary decreases in the VHS/HIMS. 
However, it is much more likely that these actually increase, as there is strong irradiation 
in the VHS/HIMS which should increase the color temperature correction, and it is not 
possible to decrease the stress on the inner boundary from the zero stress condition used 
for the HSS. Thus the data show that the VHS/HIMS is associated with a truncated cool 
disk geometry, where the disk inner radius is 1.3-2.2 times larger than that of the constant 
radius seen in the HSS which is identified with the ISCO. 



Within the Suzaku observation there is a small change in flux and spectral shape, 
and a corresponding change in the fast variability properties as determined from the 
power spectrum. The low frequency QPO increases from 4.26 ± 0.06 Hz to 5.46 ± 0.06 Hz 
as the spectrum softens, with a derived decrease in apparent disk inner radius of 
r* n = (95^g) • Cso^s km to (90^4) • Cso^s km. While this is not significant, this decrease 



in radius is consistent with the increase in QPO frequency i 



from Lens-Thirring precession of the hot fl ow interior to r\„ (llngram et al 



magneto-hydrodynamic mode of this flow (IMachida &: Matsumoto 



the QPO is forme d either 



20081 ) 



20091 ) or from a 



A slightly truncated disk in the VHS/HIMS clearly makes a smooth connection to a 
larger rad ius disk truncati on as required for LHS models which use the alternative hot flow 



solutions ( lEsin et al. 



19971 ). making a coherent picture of the outbursts of BHB in a model 
where the spectrum changes from LHS-VHS/HIMS-HSS driven by a decreasing inner disk 
radius until it reaches the ISCO. 



We would like to thank the referee Dr. Andrzej Zdziarski for his helpful comments. 
We are grateful to all the Suzaku team members. We also thank Mami Machida and Ryoji 
Matsumoto for helpful discussion on QPO and MHD. The present work is supported in 
part by Grant-in-Aid No. 19740113 from Ministry of Education, Culture, Sports, Science 
and Technology of Japan, and by Grant-in-Aid for Project Research in Shibaura Institute 
of Technology. 



-25 - 

REFERENCES 

Belloni, T., Mendez, M., van der Klis, M., et al. 1996, ApJ, 472, L107 

Belloni, T., Psaltis, D., k van der Klis, M. 2002, ApJ, 572, 392 

Belloni, T. M. 2010, Lecture Notes in Physics, Berlin Springer Verlag, 794, 53 

Cowley, A. P., Schmidtke, P. C, Hutchings, J. B., k Crampton, D. 2002, AJ, 123, 1741 

Cunningham, C. T., 1975, ApJ, 202, 788 

Done, C, Mulchaey, J. S., Mushotzky, R. F., k Arnaud, K. A. 1992, ApJ, 395, 275 

Done, C, k Kubota, A. 2006, MNRAS, 371, 1216 

Done, C, Gierlihski, M., k Kubota, A. 2007, A&A Rev., 15, 1 

Dunn, R. J. H., Fender, R. P., Kdrding, E. C, Belloni, T., k Cabanac, C. 2010, MNRAS, 
403, 61 

Esin, A. A., McClintock, J. E., k Narayan, R. 1997, ApJ, 489, 865 

Fabian, A. C, Rees, M. J., Stella, L., k White, N. E. 1989, MNRAS, 238, 729 

Fukazawa, Y., Mizuno, T., Watanabe, S., et al. 2009, PASJ, 61, 17 

Gierlihski, M., Zdziarski, A. A., Poutanen, J., et al. 1999, MNRAS, 309, 496 

Gierlihski, M., k Done, C. 2003, MNRAS, 342, 1083 

Gierlihski, M., Done, C, k Page, K. 2008, MNRAS, 388, 753 

Hynes, R. L, Steeghs, D., Casares, J., Charles, P. A., k O'Brien, K. 2003, ApJ, 583, L95 
Hynes, R. I., Steeghs, D., Casares, J., Charles, P. A., k O'Brien, K. 2004, ApJ, 609, 317 



-26 - 

Ingram, A., Done, C, & Fragile, P. C. 2009, MNRAS, 397, L101 

Kokubun, M., Makishima, K., Takahashi, T., et al. 2007, PAS J, 59, 53 

Kolehmainen, M., & Done, C. 2010, MNRAS, 406, 2206 

Kolehmainen, M., Done, C, & Diaz Trigo, M. 2011, MNRAS, 416, 311 

Kubota, A., Tanaka, Y., Makishima, K., et al. 1998, PAS J, 50, 667 

Kubota, A., Makishima, K., & Ebisawa, K., 2001, ApJ, 560, L147 

Kubota, A., & Done, C. 2004, MNRAS, 353, 980 

Kubota, A., & Makishima, K., 2004, ApJ, 601, 428 

Lasota, J.-P. 2001, Black Holes in Binaries and Galactic Nuclei, 149 

Magdziarz, P., & Zdziarski, A., 1995, MNRAS, 273, 837 

Machida, M., & Matsumoto, R. 2008, PASJ, 60, 613 

Makishima, K., Maejima, Y., Mitsuda, K., et al. 1986, ApJ, 308, 635 

Markert, T. H., Canizares, C. R., Clark, G. W., et al. 1973, ApJ, 184, L67 

McClintock, J. E., & Remillard, R. A. 2006, Compact stellar X-ray sources, 157 

Mendez, M., & van der Klis, M. 1997, ApJ, 479, 926 

Miller, J. M., Reynolds, C. S., Fabian, A. C, et al. 2008, ApJ, 679, L113 

Mitsuda, K., Inoue, H., Koyama, K., et al. 1984, PASJ, 36, 741 

Miyamoto, S., Kimura, K., Kitamoto, S., Dotani, T., Ebisawa, K. 1991, ApJ, 383, 784 
Morrison, R., & McCammon, D. 1983, ApJ, 270, 119 



-27- 

Motta, S., Belloni, T., & Homan, J. 2009, MNRAS, 400, 1603 

Muiioz-Darias, T., Casares, J., & Martmez-Pais, I. G. 2008, MNRAS, 385, 2205 

Ponti, G., Fender, R. P., Begelman, M. C, et al. 2012. larXiv: 1201 .41721 

Poutanen, J., & Svensson, R. 1996, ApJ, 470, 249 

Ross, R. R., Fabian, A. C, & Young, A. J. 1999, MNRAS, 306, 461 

Ross, R. R., & Fabian, A. C. 2005, MNRAS, 358, 211 

Ross, R. R., & Fabian, A. C. 2007, MNRAS, 381, 1697 

Rykoff, E. S., Miller, J. M., Steeghs, D., & Torres, M. A. P. 2007, ApJ, 666, 1129 

Shakura, N., Sunyaev, R., 1973, A&A, 24, 337 

Shidatsu, M., Ueda, Y., Nakahira, S., et al. 2011a, PASJ, 63, 803 

Shidatsu, M., Ueda, Y., Tazaki, F., et al. 2011b, PASJ, 63, 785 

Shimura, T., Takahara, F.,1995, ApJ, 445, 780 

Steiner, J. F., Narayan, R., McClintock, J. E., & Ebisawa, K. 2009, PASP, 121, 1279 
Steiner, J. F., McClintock, J. E., Remillard, R. A., et al. 2010, ApJ, 718, L117 
Stella, L., & Vietri, M. 1998, ApJ, 492, L59 
Svensson, R., & Zdziarski, A. A. 1994, ApJ, 436, 599 

Swank, J. H., Smith, E. A., Smith, D. M., & Markwardt, C. B. 2006, The Astronomer's 
Telegram, 944, 1 

Wijnands, R., Homan, J., & van der Klis, M. 1999, ApJ, 526, L33 



-28 - 

Yamada, S., Makishima, K., Uehara, Y., et al. 2009, ApJ, 707, L109 

Zdziarski, A. A., Johnson, W. N., & Magdziarz, P. 1996, MNRAS, 283, 193 

Zdziarski, A. A., Grove, J. E., Poutanen, J., Rao, A. R., & Vadawale, S. V. 2001, ApJ, 554, 
L45 

Zdziarski, A. A., Poutanen, J., Paciesas, W. S., & Wen, L. 2002, ApJ, 578, 357 
Zdziarski, A. A., Gierlihski, M., Mikolajewska, J., et al. 2004, MNRAS, 351, 791 
Zycki, P. T., Done, C, & Smith, D. A. 1999, MNRAS, 305, 231 
Zhang, S. N., Cui, W., & Chen, W. 1997, ApJ, 482, L155 



This manuscript was prepared with the A AS IATj^X macros v5.2. 



-29 - 



Table 1. Description of models in XSPEC 

model descriptions in XSPEC 

model 1 wabs*(diskbb+rdblur*ireflect*pl + gauss) 

model 2 a wabs*(simpl(d)*diskbb+rdblur*ireflect* simpl(c)*diskbb+gauss) 

model 3 a wabs*(simpl(d)*diskbb+rdblur*ireflect*(simpl(c)*diskbb+nthcomp)+gauss) 

model 4 a ' 6 wabs*(simpl(d)*dkbbfth(d)+rdblur*ireflect*(simpl(c)*dkbfth(d)+dkbbfth(c))+gauss) 

Note. — To use convolution models, SIMPL, ireflect, and rdblur, energies are extended in the range of 
0.1-1000 keV. 

a SlMPL(D) and SIMPl(c) represent the direct and comptonized component of the SIMPL model, respectively. 

^dkbbfth(d) and dkbbfth(c) represent the direct and comptonized component of the dkbbfth model, 
respectively. 
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Table 2. The best fit parameters for the summed spectrum 



Component Parameter model 1 model 2 model 3 model 4 



WABS 


iV H (10 21 cm" 2 ) 


6.7 ±0.2 


4.8t° i 


4.9 ±0.2 


4 9 +a i 


DISKBB 


fcT in (keV) 


0.65 ±0.02 


0.551°;°? 


0.56 ±0.02 







»*in (Csocfe km) 


56« 


1221? 


94t 6 5 


— 


NTHCOMP or DKBBFTH 


kT^ (keV) 


— 


— 


— 


U - D '-0.08 




r th 


— 


— 


2.72 ±0.04 


2 fi7 +0.02 
z - D '-0.03 




kT c (keV) 


— 


— 


r,r,+ >230 

' z -40 


41 ±6 




norm 






3.7 ±0.3 


1.2 ±0.3 




^tran (^*in) 








~ r + 0.4 

■^•°-0.2 




LIL-IIVCO / 






0.4 ± 0.3 


5S+ ' 08 




derived / th 








0.27 


PL or SIMPL 


Tpl 


2.68ir 01 


2.68 ± 0.01 


(2-1) b 


(2.1) b 




norm c 


11.1 ±0.5 










/pl 




0.33 ± 0.01 


043+ 005 


0.032 ±0.002 


IREFLECT d 


Q/2tt 


0.79 ±0.05 


n o 9 +0.03 


n oq+0.07 
0.83_ 008 


n 74+°'» 7 

u - ' 4 -0.07 




Fe abundance (solar) 


1.8 ±0.2 


1.7 ± 0.2 


1 R+ ' 4 


1 Q+0.3 
1 - !:, -0.2 


gaussian e 


E (keV) 


6.94 ±0.34 


6.96 f 


6.94 f 


6.94 f 




eqw (eV) 




o+45 
8 -8 


c+46 
°-5 


iD -16 




X 2 /dof 


190.2/186 


198.9/186 


192.0/184 


178.7/184 


photon flux g 


^dtk+Th (P^tons s-icm- 2 ) 






39.7 


43.5 


inner radius 11 


r* n ((soda km) 






117 ±5 


93tJ 9 



Note. — Errors represents 90% confidence limit of statistical errors. Energy band was extended as 0.1-1000 keV. 
a The optical depth calculated vie equation (1). 
b Fixed in the model fitting. 
c photons keV- 1 cm- 2 s- 1 at 1 keV 

d The solar abundances are assumed for heavy elements except for iron. The temperature of the reflector and 
ionization parameter are fixed at 10 6 K and 300 erg cm s -1 , respectively. The reflected components are relativistically 
blurred by RDBLUR with fixed /jrdblur = i iJ G , firdblur = w 5 Rg /3 = _ 3j and j = 50 o 

°Gaussian sigma is fixed at 0.2 keV. 
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f Value of the gaussian central energy E is limited to be 6.2—7.8. Its 90% confidence exceeds this range with 
model 2,3, and 4. 

s Unabsorbed photon flux of the disk and thermal Compton components in the range of 0.01-100 keV, where the 
reflected emission is excluded. The photon flux of diskbb plus THCOMP and that of dkbbfth are shown for model 3 
and model 4, respectively. 



h The apparent disk inner radius estimated via equation (A.l) in 



Kubota fc Makishima 



(2004) for the slab geometry. 
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Table 3. The best fit parameters of model 3 and 4 for the spectra in the first half and the 

latter half 



Component Parameter model3 model4 







first half 


latter half 


first half 


latter half 


DISKBB 


fcT in (keV) 
r- m ((50ds km) 


0.54 ±0.01 
99+J 


0.581™? 
91 ±3 






NTHCOMP or DKBBFTH fcTij?* (keV) 

r th 

kT c (keV) 
norm 
rtran (r in ) 


r, fi O+0.03 

Z - DO -0.04 

38l« 
3.5 ± 0.3 


9 sn+ 04 
> 41 

, y + 0.3 

' -0.4 


0.66 ±0.03 
2.59 ±0.03 

(41) a 

1 ,+0.3 

J ~ a -0.2 

2 6 +0 ' 4 


69+ - 02 
2.74 ±0.04 
(41) a 
1 2+ ' 1 

n q+0.4 
z '°-0.2 




derived r 
derived / th 


0.6 ±0.3 


< 0.55 


0.61 ±0.02 
0.30 


55+ - 02 
0.25 


SIMPL 


Tpl 
/PL 


(2.1) 

o rra +0 - 005 


(2.1) 
nqE+0.005 


(2.1) 
0.035 ± 0.002 


(2.1) 
0.028 ±0.002 


IREFLECT 


Q/2tt 

Fe abundance (solar) 


76+ - 09 
u - ' D -0.09 

(1.8) a 


(1.8) a 


66+ 06 
(1.9) a 


0.82 ±0.10 
(1.9) a 


gaussian 


E (keV) 
eqw (eV) 


(6.94) a 
< 49 


(6.94) a 
< 64 


(6.94) a 
< 58 


(6.94) a 

k+66 

°-S 




X 2 /dof 


185.1/184 


190.9/184 


180.5/185 


184.6/185 


photon flux 
inner radius 


^tk+Th (Photons s" 1 cm- 2 ) 
r* n ((nods km) 


38.7 
122±* 


41.1 

netl 


43.0 


44.3 


Note. — Same 


as table [2] but for the separated spectra 


fit with model 3 and 4. 







a Fixed at the best fit values for the summed spectrum in table [2] 
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Fig. 1 — Light curves of GX 339 - 4 based on the RXTE/ASM and the Swift /BAT Hard 
X-ray Transient Monitor. Panel(a) shows the 1.5-12 keV ASM count rate, while panel(b) 
shows the ASM hardness ratio (5-12 keV/3-5 keV) where data points of large statistic error 
(AHR > 0.7) are excluded. Panel(c) shows the 15-50 keV BAT count rate. A left-right arrow 
with vertical dashed lines indicates the time period during which the presented Suzaku data 
were obtained. Simultaneous observations of XMM and RXTE on 2007 February 19th is 
indicated with a downarrow. 
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Fig. 2. — The Suzaku light curves of GX 339 — 4 during the observation. Count rate of 
0.7-10 keV XISO (panel a), 10-60 keV PIN (panel b), and 50-200 keV GSO (panel c) were 
plotted against time since 2007 February 12th 05:33:31. The XISO, PIN, and GSO light 
curves are binned in 6 sec, 100 sec, and 1000 sec, respectively. The XISO events within the 
central region of 3' were included in panel (a). The PIN and GSO light curves were corrected 
for dead time, and the NXB was subtracted. A vertical dotted line divided the observation 
between the first half and the latter half. 
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Fig. 3. — The Suzaku raw spectra of GX 339 — 4, together with predictions of the best-fit 
model 1 are shown in panel (a). Residuals between the data and model 1, model 2, model 3, 
and model 4, are shown in panel (b), (c), (d), and (e), respectively. 
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Fig. 4. — vF v spectrum of GX339 — 4 from the best fit model 1, where PIN and GSO spectra 
are normaized to XIS-0 spectrum. 
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Fig. 5. — Model spectra of GX339 — 4 from the best fit model 1 (panel a), model 2 (panel b), 
model 3 (panel c), and model 4 (panel d), where the interstellar absorption is excluded. 
Emission from the underlying disk (DISKBB or DKBBFTH(d); magenta), non-thermal comp- 
ton component (PL or SIMPL(c); light blue), thermal Compton component (NTHCOMP or 
dkbbfth(c); green), and gaussian (purple) are shown solid lines. The reflected compo- 
nents of the thermal (green) and non-thermal (light blue) are shown with dashed lines. In 
panel (d), the intrinsic disk emission (magenta dash-dot line), where all the power of the 
dkbbfth component is released as optically thick disk, is plotted together with the original 
model components. 
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tran 
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Fig. 6. — The schematic geometry envisaged for the dkbbfth model. At r < r tran , sum 
of the local luminosity from thermal corona (green) and disk (magenta) ar e consistent with 



that e xpected from the standard disk at r > r tran . Details are described in 



Done fc Kubota 



(|2006|). 
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Fig. 7. — Comparison between the Suzaku averaged spectrum in the VHS/HIMS (blue) and 
XMM/RXTE spectrum in the HSS obtained on February 19 (red). vF v spectra and the 
unabsorbed model spectra by excluding the same interstellar absorption are shown in panel 
(a) and (b), respectively. A blue dash-dot line indicates the intrinsic disk emission for the 
Suzaku VHS spectrum in panel (b). The disk emission in the HSS is also indicated with a 
red dash-dot line. 
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Fig. 8. — Power spectral density plots calculated based on NXB subtracted PIN lightcurve 
in the range of 10-60 keV, where white noise level was not subtracted. Panel (a) shows the 
power spectral density based on the entire light curve, while panel (b) and (c) shows those 
of the first half and the latter half, respectively. 
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Fig. 9. — vF v spectra taken from the first half (green) and the latter half (orange) are plotted 
in the top panel based on the best fit model 4. The ratios to the best fit model 4 for the 
summed spectrum are plotted in the same color in the bottom panel. 



